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PREFACE 


The  analysis,  computer  programs  and  results  discussed  herein  are 
the  result  of  research  and  programming  performed  for 


Rome  Air  Development  Center 

Electromagnetic  Sciences  Division 

L.  G.  Hanscom  Air  Force  Base,  Massachusetts  01731 


mm 


1 .  INTRODUCTION 


During  the  last  30  months,  our  technical  work  consisted  of  providing 
RADC  scientists  with  mathematical  and  programming  support  in  the  following 
areas:  magnetostatic  wave  transducer  analysis;  antenna  analysis  and  design 
studies;  conversion,  implementation,  and  execution  of  ionospheric  propaga¬ 
tion  programs;  radar  data  analysis;  electromagnetic  scattering  from  rough 
terrains,  etc.  This  work  resulted  in  the  development  of  mathematical  models 
and  computer  programs.  When  desired,  we  functioned  in  a  collaborative  manner 
to  assist  with  the  physical  formulation  of  problems  of  interest.  This  report 
presents  a  representative  selection  of  our  efforts. 


2.  MSSW  Transducer  Analysis 


2.1  Introduction 

The  development  of  mathematical  models  for  describing  the  response  fea¬ 
tures  of  magnetostatic  surface  wave  (MSSW)  devices  has  been  of  considerable 
interest  to  the  Air  Force.  Despite  previous  successes  at  predicting  the 
behavior  of  transducers  lifted  off  a  YIG  substrate,  the  models  have  proved 
to  be  inadequate  for  describing,  with  sufficient  accuracy,  the  response  of 
multi-element  transducer  arrays  very  near  or  on  a  magnetically  biased  YIG 
slab.  Strong  coupling  effects  appear  to  become  important  in  these  situa¬ 
tions.  Our  efforts  during  this  contract  were  aimed  at  developing  models 
for  predicting  transducer  characteristics  under  strong  coupling  conditions. 

In  keeping  with  the  initiator's  objectives,  a  sequence  of  investiga¬ 
tions  was  undertaken  to  develop  current  density  models  for  multi-element 
arrays  on  or  near  a  YIG  slab.  The  first  study  entailed  determining  current 
density  over  mutually  coupled  conducting  strips  near  a  ground  plane;  no 
YIG  material  being  present  to  influence  the  distributions.  The  Fourier 
transforms  of  the  resulting  current  profiles  were  incorporated  into  the 
Weinberg  programs  [Ref.  1]  for ‘predicting,  among  other  things,  transducer 
insertion  loss.  The  details  of  the  analysis  as  well  as  the  results  have 
been  previously  described  in  Ref.  2.  The  current  distributions  did  not 
lead  to  some  of  the  insertion  loss  features  that  appear  experimentally. 

Even  when  replacing  the  ground  plane  with  a  highly  permeable  (u  -*•  00 ) 
magnetic  half-space,  the  insertion  loss  predictions  did  not  improve. 

From  the  investigations  it  was  concluded  that  mutual  coupling  among  the 
strips  and  ground  plane  (or  magnetic  half-space)  does  not  by  itself  appear 
to  be  the  key  mechanism  for  causing  some  of  the  dominant  response  features 
observed  in  the  surface  wave  passband. 


The  results  noted  above  led  us  to  several  Green’s  function  analyses. 
By  posing  (and  solving)  a  sufficiently  well  modeled  boundary-value  problem 
describing  the  field  about  a  time-harmonic  line  current  source  near  a  TIG 
slab,  in  principle  one  can  obtain  the  current  density  on  a  planar  array 
of  infinitely  long,  thin,  perfectly  conducting  strips  on  or  near  a  ferrite 
slab.  The  current  density  on  each  strip  can  be  thought  of  as  arising  from 


a  continuum  of  line  currents.  Obtaining  the  Green's  function,  or  total  in¬ 
fluence  function,  of  the  magnetically  biased  slab  on  the  vector  potential 
field  about  the  line  source  became  the  focus  of  most  of  the  work  to  be  re¬ 
ported  in  this  section. 

2.2  Transducer  Green's  Function  Analysis  -  Introduction 

The  geometry  of  MSSU  transducers  is  often  such  that  source  and  field 
quantities  can  be  assumed  to  vary  much  more  strongly  with  x  and  y  than 
z.  A  two-dimensional  field  theory  approximation  was,  therefore,  invoked 
for  determining  the  Green's  function  of  an  array  between  a  ground  plane  and 
YIG  slab.  (See  Fig.  2.1.)  The  strips  were  also  assumed  to  be  sufficiently 
close  to  the  YIG  and  each  other  so  that 

o  <  <<  ; 

is  the  electromagnetic  wave  number  in  the  medium  (GGG)  above  and  below 
the  YIG,  and  R  is  the  distance  from  any  line  current  element  on  a  strip 
to  any  field  point  of  interest.  Field  points  on  nearby  strips  or  on  the 
nearest  regions  of  the  YIG  slab  surfaces  are  most  important  regarding  mutual 
coupling  and  generation  of  MSSWs.  It  was  assumed,  therefore,  that  all  field 
points  required  for  a  current  density  analysis  were  located  in  the  quasi¬ 
static  zone  of  the  rf  field  of  the  line  currents.  At  frequencies  between 
2.5  -  3.5  Ghz,  this  assumption  appears  to  be  reasonable  for  typical  MSSU 
arrays  composed  of  a  few  strip  elements. 

2.3  Derivation  of  Two-Dimensional  Field  Equations 

The  two-dimensional  approximation  discussed  above  leads  to  a  consider¬ 
able  simplification  of  the  analysis.  Maxwell's  equations  must  be  satisfied 


in  the  3  regions  as  must  the  boundary  conditions  at  the  surfaces  y 
y  *  0,  and  y  *  -h.  With 

_,£>  ,  0>  .O 

s  -/»•  H  =  V  X  A  ; 

fcf .  y  x  f  > 

tP  *  i 

V  •  T2>^  -  5  •  =  V  *  o  ^ 

<7-7^  --  ; 


(based  on  the  linearized  magnetization  equation) 


for  the  surface  wave  case  in  which  y  ■  y  .  Maxwell’s  equations  in  the 

11  22 

3  regions  take  the  following  form: 

C  (?*  A®)  s/t.  J  ; 

V  *■  (/£~‘  V  <  A®  )  =  fJo  ; 

A*)  ; 

1<  E®  -  -  Z(v*f). 


( 9 )  yields 


and  ( 8 )  becomes 


/  \  V7a  T 

V  (l *  A  )  -  V  A  =  a,  0  +  cyv. 

Inserting  the  expression  for  into  the  above  gives 

Y(V-A®)-^lA 4 


Applying  the  Lorentz  condition: 


V-  A  +*  —  & 


leaves 


Similarly,  for  region  (3): 


r  a 


when 


A<2> 

V*  A  4-  ^JA0  ^ 


N^VSV.V^V^-V.>:.:' 


‘  «*  y*  \*  * 


' -  V»»>  j 


For  region  (2),  since 


v  x.  [fa'  v  x  A®)  =  (v- 1 f)  -  fa fa 


we  get 


fa  i  fa  fa)  -  fa  fa  3.)  fa- -■ /*€^Sfa§ 3 


The  desired  guage  requires  the  condition 

fa 7  fa  fa)  =o} 


leaving 


Since  J  and  H  are  aligned  parallel  to  the  z-axis  , 


Assuming 


A,  ■  A  e- 

and 

T*--  3 e* 

the  exact,  linearized  wave  equations  governing  the  vector  potential  fields  in 
the  3  regions  reduce  to 

A  =  -  ft.  J  ; 

V  )  A®  ^  4s-  o  ; 

A®  +  /tc<AlaL  A&  =  o. 

Expanding  the  wave  equations  in  Cartesian  components  leads  to  the  following 
equations  for  a  line  current  at  x  ■  0,  y  -  d: 


US 


A  =  -y/ol  £(*0  St y-J) 


4<P  .© 

+  4 


Af^  A~=<3  j 


© 


A«  4  A®.  + 


yv 


/‘o€y  >4^  -O  . 


As  noted  earlier,  the  current  on  any  strip  can  be  expressed  as  a  continu¬ 
ous  density  of  line  currents.  It  follows  that  the  total  vector  potential  for 
a  strip  can  be  regarded  as  a  superposition  of  the  vector  potentials  caused  by 
a  weighted  density  of  line  currents.  The- vector  potential  for  a  single  line 
current,  i.e.,  the  Green's  function  of  the  system,  can  be  found  by  solving 
the  system  of  wave  equations  subject  to  the  boundary  conditions.  The  Green's 
function,  G,  is  proportional  to  the  vector  potential  and  contains  all  of  the 
influences  that  affect  the  current  densities.  The  boundary  conditions  that 
A  (or  G)  must  satisfy  will  now  be  derived. 

Since 


it  follows  that 


with  m  ■  1,3  for  regions  1  and  3.  For  the  slab  region 


The  boundary  conditions  at  the  ground  plane  and  YIG  surfaces  may  be  summar 
ized  as  follows: 

so  erf  . 

.g©  =  ^  j 

A  ^ 

,  ,G>  ,  .<& 

H*  *  H*  A+  V=0  y 

•e£  -  ..j 

Uf  =  Mf  *4  y  s  -  K  . 

Using  the  relationships  among  IJ,  H  and  A,  the  boundary  conditions  may  be 
written  in  terms  of  the  A  ®  alone: 


V.' 

V-1 

.1 

1 


Y*  ? 


txi'  y=o 


a© 

A. 


A®  --  A® 


(45) 
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.  c  -VJ 


The  propagation  constants  for  the  3  regions  are  given  by  equations  (29-31): 


C  ■=  jUt, 


(49) 


(50) 


and 


£  =  yMo  G3  lO. 


(51) 


At  frequencies  within  the  surface  wave  passband,  k2  is  negative  and  is 
denoted  by  k2  =  — k2.  In  the  quasi-magnetostatic  approximation  the  wave 


equations  < 

are 

approximated  by 

A® 

A  ® 

Z  - yi/0  jT  S'fXJ 

+ 

A  yy 

>41 

+ 

A® 

r\  yy 

z  C 

j 

A® 

A® 

A 

z  o  . 

(52) 


(53) 

(54) 


This  is  equivalent  to  neglecting  the  displacement  current,  i.e.,  radiation 
term,  in  the  curl  H  equation  for  each  region. 


2.4  Extended  MSW  Analysis 

We  found  that  it  was  important  to  retain  the  term  involving  (k2  A  © ) 

s 

for  the  YIG  region.  That  term  has  very  little  effect  on  the  MSW's  wave- 
n umber  over  most  of  the  surface  wave  band;  however,  it  introduces  a  long 
period  surface  wave  in  addition  to  the  MSSW  and  other  Important  effects. 
The  extended  quasi-magnetostatic  surface  wave  equations  may  thus  be  sum- 


marized 

as  follows: 

.<Z> 

Cb 

AxX 

+  Ayy 

A& 

A<S> 

+  Ayy 

-  fcr  A® 

A® 

A® 

*  Ayy 

<3  . 

Upon  solution,  B  and  E  follow  from  (2,14)  and  similar  formulas  for  the  other 
regions:  since  V<£  ■  ^  ■  0  in  a  pure,  two-dimensional  formulation  where 

there  is  no  z-variatlon,  <p  is  not  required  in  the  expression  for  _E. 

The  solutions  for  each  region  that  satisfy  the  extended  MSW  equations 
may  be  written  as: 


a  £  2  A  L*v+t"J>‘3  +  e  U[s,e'‘ 

-oO 


-ifelT 


f Umx/4-") 


ao 

t  J,  y  -y <4^+41  1 

=  jdt  e  [C,e  +  Cte  J; 
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a 

a 


- 


(f+k)l*-l 

Jk  €.  X),  <S  # 


For  region  the  direct  field  due  to  the  line  source  alone  has  been 

separated  out  of  the  general  solution,  i.e.,  the  unbounded  Green's  function 
term 


-Qih  JL  L*’+ 


=  <-2^L  JL  [**+  • 

All  potentials  are  then  made  dimensionless  by  dividing  out  the  factor  (i^lMir) 
for  convenience.  The  normalized  vector  potential  or  Green's  function  for 
each  region  is  to  be  found  for  a  line  source  at  (o,d).  The  |k|  -  factors 
in  the  expressions  for  A®  and  A®  are  rigorously  correct  only  for  pure 
static  fields.  The  implications  for  a  quasistatic  analysis  become  apparent 
shortly.  Also,  while  the  Green's  function  for  a  quasistatic ‘line  source 
varies  as 

J. L  k,~R  UjY-J  ) 

=  JL  4.  +  , 

the  term  £n  kj  may  be  dropped  now  (or  later)  because  all  the  boundary  con¬ 
ditions  involve  spatial  derivatives  in  x  and  y.  The  vector  potentials 
are  thus  arbitrary  to  within  a  constant  level  of  potential.  Bj,  C^, 

Cj,  and  are  all,  as  yet,  unknown  functions  of  k.  Inserting  the  solu¬ 
tions  into  the  boundary  conditions,  and  solving  for  the  unknown  spectral 
amplitudes,  we  arrived  at  the  following  expression  for  a®  in  the 


important  limiting  case  for 


where 

aXt)  =  ctx  -  uicr+J), 

=  iw[>- 

a  if  H 

and 

©CM  =  [lkl%  k'CeV»v)+  kll'3 

+  zEIMVFTv  cM,  <t\TSr+3T 


One  of  the  findings  of  the  extended  MSW  formulation  was  that  £(k) 
zeros  at 


k  =  ±  ke 


and 


possesses 


The  poles  of  the  Integrands  at  t  k0  are  very  near  the  origin  and  lead  to 
very  long  period  surface  waves.  Earlier  MSSW  analyses  do  not  reveal  such 
poles  because  of  their  simpler  characteristic  equation;  namely. 


©  ( *0  =  0  =  C1*  V)  +  2 r  CrtL  C  I klh). 

ks-+  o 


The  poles  near  the  origin  are  located  to  an  excellent  approximation  at 

ko  &  -  i  1  1 1  K  • 

Gerson  and  Nadan  [Ref.  3  ]  report  other  "dynamic"  modes  in  their  plane  wave 
analysis  involving  the  exact  set  of  Maxwell's  equations.  A  comparison  be¬ 
tween  their  results  and  ours  in  this  regard  is  by  necessity  incomplete  as 
of  this  writing. 

The  problem  that  remains  for  determining  the  Green's  function  is  the 
evaluation  of  the  two  Fourier  integrals  shown  in  (  63  )  •  They  contain  the 
effect  of  the  YIG  on  the  total  vector  potential  and  are  denoted  by 

CD 

§oo 

oC 

(S  -  • 

Ifcf  S(M 

Since  b  k*  *  ^(o) ,  and  taking  into  account  that  with  an  exp(-it)t)  varia¬ 
tion  fields  propagating  in  the  ±  x  directions  correspond  to  exp(±ikx). 
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as  follows: 


we  can  express 
00 


aL  =  iJt 


±t*fex.  .)KUy-W) 


o 

co 


+  Jit  el 


kx  - 

£. 


l)t|  £(*0 


(72) 


*  a± 

is  to  be  evaluated  for  a  field  point  at  +x  and  source  point  at  x=0; 
CL_  for  the  field  point  at  -x  and  source  at  x*0.  for  a  low  loss  YIG 
slab,  p  (k)  possesses  zeros  at  k  =  ±  k0,  ikjjggy*  none  of  which  are  at  the 
origin;  all  of  which  are  very  near  the  real  k  axis. 

The  2^  integral  in  particular  requires  special  treatment  as  will  now 
be  discussed.  Note  that  it  is  equal  to 


eO 

f  £»**  _  k(yw)  r 

e  c  | 

■  b<*> 

*  lfe.1  v 

_  BtM 

eO 

"] 


JL  I 

a  o  J 


it’tx  -fc.  (.Y+J) 

dt  e  e 

!  | 


(73) 


The  1 —  part  is  non-singular  at  k*0,  and  the  2 —  may  be  related  to  the 
Exponential  integral  of  complex  argument.  Let  k  ■  ev  .  It  follows  that 
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.A 


e  € 

=  JL  L,  cel” 

C-»o  t  If 

SJL  £.le<y+J 
£  *»« 


See  Ref.  4 


For  e-*o. 


with  Yg  “  Euler's  constant.  Since  spatially  uniform  terms,  however  large  or 
small,  contribute  nothing  of  physical  consequence  to  the  vector  potential, 
such  terms  may  be  discarded.  Collecting  physically  meaningful  terms  leads 


Qt=rfc  e 


DM 

u 


©CM 

.tO+j-.-ixn  _ 

e  f  °‘°> 


l! 


-  I 


-  [V+  ^ 


For  an  analysis  based  on  the  exact  equations  (32-34  )  and  =  e3  ,  |k| 
would  be  replaced  by  /k*  -  k|  .  In  the  quasistatic  approximation,  one  must 
think  in  terms  of  propagating  fields  with 


—  yL.  -  H-i'-j 


k,-*o 


C  *o 


hence  the  need  for  the  limiting  process  in  (  73  )  .  For  a  pure  static  analysis, 
kj  *  0  identically.  No  resolution  into  separate  propagating  fields  in  either 


the  +x  or  -x  direction  is  to  be  made  in  this  case.  One  encounters, 
rigorously,  an  integral  of  the  form 


jJc-X  . 

e. 

iM 


that  holds  for  all  x  and  y  and  whose  "value",  in  the  sense  of  generalized 
functions,  is  given  by 

—  JLi  0/+jy^  4 

See  Lighthill  [Ref.  5  ].  Assembling  terms  and  making  the  following  change 
of  variables,  the  integrals  comprising  are  then  ready  for  evalua¬ 

tion.  Let 

feM  -  kns sw, 

(J  *  k.  ^  J 

4  =  L/+c/ 

=  kM  [y-w-t  i <.» 3 /C «k« k ^ 

s  [  |3  -  i  1*0  /  • 

The  total  vector  potential  becomes 

<3°  «  JL  L*l+  +  at 


where 

a. 


*  i  l  ©(y>  u  J 


with 


£)(y)=  Q(°)  +  v'il+e+Z')  tzlu  iu'+$  Cf&t'Iu^+oj’ ) 


(92) 


and 


B(.)  =  fc'tV  r 


(93) 


2.5  Evaluation  of  CL  ± 

A  direct  numerical  evaluation  of  CL  +  was  accomplished  using  the  methods 
of  Real  Axis  Integration  [Refs.  6,7].  For  a  nearly  lossless  medium  there  are 
2  poles  along  the  path  of  integration  which  is  infinitesimally  near  the  real 
O-axis.  The  pole  locations  are  at 


u  =  a  e  4  )  l  \  V* 

(J  r  S 


(94) 


(95) 


where  U*kh,  0o  -  k0h,  and  0  *  k  h.  To  evaluate  the  one-sided  Fourier  inte- 

s  s 

grals  in  the  sense  of  Cauchy  Principal-Value,  the  semi-infinite  range  (o  ,  °°) 
is  partitioned  into  the  ranges  (0,2U.,)  and  (2Uvf,00).  The  two  poles  are 
located  in  the  first  interval  and  satisfy  the  relationship 


O  <  U.  «  . 

The  problem  reduced  to  the  evaluation  of 

•O 


eo  &<•> 

T  .so .  f  &«>  Ilf 

Jju  £  Vtui:  jjoe  |_£,0>  +  U  j 


o 


00 


=  -fro) 


■20* 


(96) 


(97) 


(98) 
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The  singularities  at  Uq  and  were  analytically  removed  from  the 

integrand,  and  their  contributions  added  back  into  the  final  answer  separately, 
viz. , 


a 

O  O 

a 

f  _yu  [~  C>  _  n 

+  J U-t;o  u— j 


(99) 


The  2 —  integral,  although  singular,  is  easily  evaluated;  the  1 —  integral  is 
now  non-singular.  To  determine  CQ  and  CM>  express  f(u)  as 


$ 


CO) 


ncv) 

8 1»)  ’ 


where  n(u)  is  a  numerator  function. 
Then,  as  u+u 


Co 

?  (O)  =  U 


(100) 


(101) 


or 

c =  JL 


It  follows  that 


A 


u-«-  t> 


*  ©(*>•)  (jLzXslL  ®<Uo)  4 - 


n  <*•>») 
&'<v '•> 


since  0(Uq)  *  0  at  the  pole.  Similarly, 


_ 

c-m  — ;p  .  • 

A  32-point  Gauss  quadrature  evaluation  of  the  nonsingular  integral 

_ 

[?,».-  fe.  -  ] 


compared  well  with  64-point  Gauss  quadrature  evaluation  in  the  sense  that 

discrepancies  between  the  two  appeared  at  the  10^  significant  digit.  It 

proved  crucial  to  evaluate  C  and  Cu  using  analytical  expressions  for 

o  rl 

jB'(U).  Finite  difference  approximations  gave  poor  results.  The  evaluation 
of  the  singular  integrals  must  now  be  considered. 


One  can  express  the  1 —  of  the  singular  integrals  as  follows: 

T  r  °  -*>•  c.  t,  c. 

\JO&  -JTo.  *  3Joe  +]dtJe  <"«'• 

A  O 


[  ,  -*>  C. 

]dt,«  - 


Taking  into  account  that  |yU  |«  1  for  any  conceivable  transducer,  and 
U  -  U  0 

letting  W  *  — - -  ,  one  is  led  to  the  excellent  approximation  that 


au0 


5du  e 


-VU  Co 


+-- 


—  Vo 


For  the  2 —  integral 


\ao  «<l£. 

J  U  — 


1+1*0 


(aoM-  u.)/ua 

=  C.  (dw  c.TO‘C 

J  w 


=  c.  e^°\E,  1-tV'l  -  £' [*(*4-**)] j 


(110) 


(111) 


Collecting  terms  one  finally  arrives  at 

.  f 


i  -iv 

lu  e  — 


E,l?U0  "  E«y(2t;A1-y,,^]je.  (112) 


The  other  singular  Integral  was  cast  into  the  form 


Om  f  VcuVY 


(113) 


0  U-0%  -* 

which  was  easily  and  accurately  evaluated  using  an  even-point  Gauss  quadrature. 
This  completes  the  evaluation  of  the  integral 

(  .  -Vo  r 


do  e"w^u,. 


(114) 


2.7  Integration  Over  the  Semi-infinite  Interval  (2U-,,°») 

For  a  line  current  residing  on  the  YIG  film,  the  integrand  can  fall  off 
as  slowly  as  1/U  for  large  U.  We  therefore  "subtracted  off"  the  asymptotic 
behavior  of  the  integrand  and  then  "added"  its  effect  brfck  into  the  final 


answer  as  shown: 


“do  e^-U 


-  Jjo  e  *>0  -  £-  Uu"v-1 

2.(Jm  . 

oO 

to 

f  / t; 

=  Jdu  eVU[vuJ'%JWJ]  +  e  %>'u) 

auM  iuh 


-  E,  lio~n  3  . 


^^(U)  is  an  asymptotic  form  of  (U)  and  its  use  as  indicated  above  speeds 
the  convergence  of  the  integral.  Since 


s  *  ~b~ 


=  - 

©Co* 


it  follows,  after  some  algebra,  that 

?«£u>  =  4  S«-»o*'3  S, 


- 1 
60 


£«  =  (i-  1 1 0  it  1  ^ 


and 


8*  =  C  •+*'+ 


It  can  be  seen  that  ^^(U)  is  expandable  in  partial  fractions  and  that  one 
is  led  to 


Jdo  =  F,  E,  1^3 


>3 


+  <2  £, 


where 


F.  = 

fv  = 


-  Sf»>  /  £o 

i[^6»  +  J 


and 


F,  =  fO 


zr=n 


2.8  Half-Period  Integration 

The  only  integral  remaining  to  be  evaluated  is 

aO 

X  *  \ju  ^ 


ao- 


Let 


Of°>  =  -  %»cu> 


(122) 


(127) 


(128) 


and  replace  y  with  (8  -  i|x|)/U^  in  the  integral. 

This  leads  to 

o #  — 

J.  Lj)  &  QLu'J’O*  1 

'*•  i 

The  procedure  adopted  for  its  evaluation  was  to  integrate  over  N  half-periods 
of  the  sine  and  cosine  functions  integrals  separately.  4  sets  of  partial  sums 


were  thus  obtained 

oO 


left  half-period  integrations 
right  half-period  integrations 


left  half-period  integrations 
right  half-period  integrations 


Each  sequence  of  partial  sums  was  then  extrapolated  to  N  +  00  using  a  Shank's 
transformation  [Ref.  9  ].  A  modified  Shank's  transformation,  was  actually 
implemented  [Ref.  10]  but  it  did  not  do  better  than  the  simpler  version. 


<f  5  ©  +  Cn-i)-TT 

O*  0  <  TT* 


*  I*  -  -  •  j  ^ 


it  follows  that 


ns  *  © 


t  \_G+ 


showing  that  each  term  in  the  partial  sum  involves  a  half-period.  For  the 
actual  computations  over  each  half  cycle,  it  was  useful  to  let  0  =  (t  +  1) 


.V'-T-WT 
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2.9  Results  and  Discussion 

The  quantity^®,  the  normalized  Green's  function,  was  programmed  and 
put  through  a  series  of  tests.  As  Fig.  2.3-2.21  shows,  the  Green's  function 
has  real  and  imaginary  parts  as  depicted  by  solid  and  dashed  lines  respectively. 

The  real  and  imaginary  components  are  off  by  a  physically  irrelevant  complex 
constant  in  the  figure.  The  amplitude  in  the  -x  direction  is  larger  than 

a 

in  the  +x  direction,  as  would  be  expected  based  on  the  direction  of  x  n. 

See  Fig.  2-1  and  [Ref.  8].  Furthermore,  the  ±x  amplitudes  are  not  uniform. 

From  the  tests  we  were  able  to  perform,  this  does  not  appear  to  be  numerical 
error:  there  seems  to  be  a  small,  spatially  dependent  transient  near  the 
line  source  at  the  origin.  At  large  distances  from  the  source,  harmonic  surface 
waves  dominate;  very  near  the  source,  as  seen  in  Fig.  2.2,  there  is  a  logarithmic 
variation  in  the  field,  as  well  as  discontinuity  at  the  source;  and  between  the 
two  extremes  a  region  of  transition. 

For  distances  of  150  microns  or  less  from  the  source,  and  using  the 
same  parameters  noted  earlier,  the  integral  I  in  (127)  was  found  to  be 
negligible.  At  considerably  larger  distances  the  half-period  integrations 
and  extrapolations  become  important:  the  algorithm  showed  convergence  to 
9  or  10  significant  digits  at  distances  of  about  A„  or  less;  whereas,  at 

n 

3  Ajj  6-8  significant  digits  were  obtained  for  I  at  several  frequencies. 

For  large  distances,  pure  harmonic  surface  waves  of  the  form 

StG  (138) 

strongly  dominate  the  Green's  function. 

The  usefulness  of  the  Green's  function  for  calculating  the  surface  cur¬ 
rent  distributions  on  one  or  more  strips  depends  to  a  large  measure  on  the 


time,  i.e.,  expense,  of  its  computation.  For  a  given  frequency  and  a  given 
distance  (positive  or  negative)  of  the  field  point  from  the  line  source,  it 
took  less  than  0.15  seconds  to  evaluate  the  Green's  function.  Additional 
programming  refinements  and/or  a  loosening  of  some  of  the  accuracy  tolerances 
would  reduce  the  computation  time  even  further;  however,  the  program  appears 
to  be  a  useful  tool  in  its  present  form  especially  when  used  in  conjunction 
with  a  cubic  spline  interpolation  code.  The  next  stage  of  work  is  to  imple¬ 
ment  the  Green's  function  for  obtaining  the  current  density  on  a  strip  re¬ 
siding  on  the  YIG.  This  involves  inverting  an  integral  equation,  and  is  in 
progress. 


Fig.  2.2  Green's  Function  at  f  =  2.8  Ghz 


2.16 


Fig.  2.21 
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MODESEARCH  and  FASTMC  are  a  tandem  set  of  computer  programs  written 
in  the  UNIVAC  version  of  ASCII  FORTRAN.  ARCON  successfully  converted  these 
programs  for  use  on  the  CDC  6600  computer.  A  full  description  of  the  pro¬ 
grams,  including  the  theoretical  background,  may  be  found  in  References  [1,2]. 

The  MODESEARCH  program  is  guided  by  control  cards  which  allow  the  user 
to  specify  ionospheric  profiles  and  ionospheric  collison  frequency  profiles. 
The  FORTRAN  NAMELIST  definition  is  used  to  enter  other  specifications,  e.g.. 
the  geomagnetic  field,  the  ground  conditions,  the  distance  of  this  slab 
from  the  transmitter.  [See  Appendix  I] 

Various  output  options  are  available.  The  FASTMC-compatible  output  is 

a  card  image  file,  in  which  the  first  card  identifies  the  slab  conditions. 

It  is  followed  by  the  mode  constants,  two  cards  per  mode,  ending  with  a 

blank  card.  Appendix  I  illustrates  a  typical  stream  of  control  code,  input 

data  and  output  data. 

• 

FASTMC  implements  a  mode-summing  technique  which  calculates  field 
strength  at  user-specified  transmitter  and  receiver  altitudes.  [See  Ap¬ 
pendix  II  for  control  cards  and  sample  input/output  data].  Output  options 
include  the  vertical  and  horizontal  components  of  the  amplitude  and  the  cor¬ 
responding  distance  from  the  transmitter,  the  results  of  which  may  be  printed 
or  plotted.  In  this  instance,  plots  were  obtained  by  using  the  TEKTRONIX 
Graphics  Terminal  (off-line  mode)  and  its  Hard  Copy  unit.  [See  Appendix  III 
for  examples  of  plots] . 
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APPENDIX  I 


<TOP  OF  FILE> 

HOLST,T300,  CM  20000.  3334  HOLST 

ATTACH, XX, NODE SEARCH, I  D  =  HOLST. 

FTN5( I=XX,L=0) 

REQUEST, TAPE1 ,*PF. 

L60. 

CATALOG,  TAPE1,«ODESEARCHOUTPUT,ID=HOLST,HR=1. 

RETURN, XX, LGO,T APE 1 . 

EOR 

NAME 

IDATUn 

AZIH=283.0,  C001P-21 .0,  HAGFLD=5.44E-5, 

EPSR=1 0.0,  SIGMA=1 0.0E-3, 

BETA=0.5,  SCLHTS-6.0,  ENHIN=0.1,  HPR1ME=87.0, 

FREQ=23. 0,  RHO=0.0,  REFLHT=50.0, 

RANGER-70.0,87.0,  RANGEI=0.0,-1 .0, 

TEND 

QUIT 

< BOTTOM  OF  F1LE> 


<TOP  OF  F1LE> 

R  0.000  F  23.0000  A  283.000  C  21.000  If  .544E-04  S  1.000E-02  E  10.0  T  87.0 

1  89.91283  -6. 502712-3. 33252104E-05-2.37283889E-04-1  .29396032E-09-1  .55444857E-1 0 

2  89.91283  -6. 502712-4. 1 3669031 E-07-3.78898638E-07  9.90072854E-01  6.71 432292E-02 

1  89.74584  -5.505621  5.50719615E-04-4.78805072E-03-6.42842801E-10-5.42580297E-11 

2  89.74584  -5.505621  1 .16720717E-06  1 . 3322761 2E-06  9.88844736E-01  6.61828827E-02 

1  87.88375  -.913771-5.32683379E-03-2.30771489E-02-4.32628737E-09-3.35977677E-09 

2  87.88375  -. 913771-6. 48275418E-06-9.45516889E-06  9. 85049639E-01  6.38126194E-02 

1  85.91200  -.370342  4.24848636E-03-8.75920738E-03-2.23193224E-08  2.3831 9225E-09 

2  85.91200  -.370342  8.94973443E-06  1 . 1 901 1255E-05  9.83474238E-01  6.30287850E-02 

1  81.32110  -.445981-5.66746234E-03-1.91768894E-02-2.37875909E-08-1 .73385229E-08 

2  81.32110  -. 445981-1. 47369354E-05-1 .961 20402E-05  9 . 76634392E-01  5.82415272E-02 

1  80.10598  -.303952  5.81933423E-03-9.44696655E-03-6.3351 1608E-08  1 .605673B8E-08 

2  80.10598  -.303952  1 .67737623E-05  2. 1 4987653E-05  9.73579686E-01  5.6729441  IE-02 

1  76.27959  -. 359672-7. 38935438E-03-1 . 01 423447E-02-1 .01 074771 E-07-5.39628050E-08 

2  76.27959  -. 359672-2. 97935702E-05-2. 461 64902E-05  9.61 396854E-01  4.85706901E-02 

1  75.23709  -.451541  7.98421 382E-03-1 .59733251E-02-7.39162419E-08  5 .26071 924E-08 

2  75.23709  -.451541  3. 1 8594953E-05  2. 59897420E-05  9.5731 6804E-01  4 .52404947E-02 

1  71 .59467  -.281 662-5. 24635889E-03-3.32464239E-03-2. 501 951 09E-07-6.83767739E-08 

2  71.59467  -. 281662-3. S1979008E-05-1 .60890849E-05  9.371 801 59E-01  3.21704217E-02 

< BO T T 001  OF  FILE> 
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APPENDIX  II 


EDIT. 

Pi 

<TOP  OF  F1LE> 

HOLST, CM1«0fi00f T160. 

ATTACH, XX, FASTMC,IJ=HOLST. 
FTN5(I=XX,L=«) 

ATTACH, TEK,TEKOFFLINE. 

LIBRART(TEK) 

ATTACH,  TAPES,  FASTJICBATA,  1D=H0LST. 
Ft£PU£ST,T4*140, *PF. 

LOO. 

CATALOG,  T4B140,FASTMCPL0T,ID=H0LST. 
RETURN, TAPES. 

RETURN, XX.LGO. 

<BOTTOH  OF  FILE> 


3334  HOLST 


APPENDIX  II  (continued) 


<T0P  OF  FIL£> 

NAME 

JDATUH 

NPRINTS1 ,  1C0J1P=1,  NRCURV=2, 

TALT=0.0,  RALT=f .0, 

A«P«IN=-.7£*02,  SIZEY=7.0, 

1END 

DATA 

FASTflC  TEST  RUN 

R  0.000  F  23.0000  A  283.000  C  21.000  it  .544E 

1  89.91283  -6. 502712-3. 33252104E-05-2.37283889E 

2  89.91283  -A. 50271 2-4. 1 34698 31 E-07-3.78898638E 

1  89.74584  -5.505621  5.50719615E-04-4.78805072E 

2  89.74584  -5.505421  1 . I4720717E-06  I.33227612E 

1  87.88375  -. 913771-5. 32683379E-03-2.30771489E 

2  87.88375  -. 913771-6. 4827541 8E-06-9. 4551 6889E 

1  85.91200  -.370342  4.24848636E-03-8.75920738E 

2  85.91200  -.370342  8.94973443E-06  1.19011255E 

7  81.32110  -. 445981-5. 68746234E-03-1 .91768894E 
2  81.32110  -. 445981-1. 47369354E-05-1 .96120402E- 

1  80.10598  -.303952  5.81 933423E-03-9.44696655E- 

2  80.10598  -.303952  1 . 47737423E-05  2.14987653E- 

1  74.27959  -. 359672-7. 38935438E-03-1 .01423447E- 

2  76.27959  -. 359672-2. 97935702E-05-2.46164902E- 

1  75.23709  -.451541  7.98421382E-03-1 .59733251E- 

2  75.23709  -.451541  3. 18594953E-05  2.59B97420E- 

1  71.59467  -. 281662-5. 24635889E-03-3.32464239E- 

2  71.59467  -. 281662-3. 81979008E-05-1 .60890849E- 


1  80.10598 

2  80.10598 
1  74.27959 


-04  S  1.000E-02 
-04-1 .29396032E 
-07  9.90072854E 
-03-6.42842801E 
-06  9. 88844736E 
-02-4.32628737E 
-06  9.85049639E 
-03-2.231 93224E 
-05  9. 83474238E 
-02-2.37875909E 
-05  9.76634392E 
-03-6.3351 1608E 
-05  9.73579686E- 
-02-1.01074771E 
-05  9.61396854E- 
-02-7.391 6241 9E- 
-05  9.57316S04E- 
-83-2.50195109E- 
•05  9.37180159E- 


E  10.0  T  87.0 
-09-1. 55444 B57E-10 
-01  6.71432292E-02 
-10-5.42580297E-11 
-01  4.61 828827E-02 
-  09-3. 35977 677E -09 
-01  6.38126194E-02 
-08  2.3831 9225E-09 
-01  6.30287850E-02 
-08-1 .73385229E-08 
-01  5.82415272E-02 
•08  1 . 60567308E-08 
-01  5.6729441  IE-02 
07-5.39626B50E-08 
-01  4.85706901E-02 
08  5.26071924E-08 
01  4.52404947E-02 
■07-6. 8 376773 9E -08 
01  3.21704217E-02 


3.000  F  23.0000  A  283.000  C  21.000  N  .544E-04  S  1.000E-02  E  10.0  T  87.0 


R  7.000  F  23.0000  A  283.000  C  21.000  M  .544E-04  S  1.000E-02  E  10.0  T  87.0 


R  40. 

NAME 

IDATUfl 

1C0HP»2, 

SEND 

START 

<8077011  OF  FILE> 


L%*.« 
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.  OTH  Radar  Data  Analysis 
4.1  Introduction 

The  Over-the-Horizon  Experimental  Radar  System  (ERS)  operated  from 
June  1980  to  June  1981.  Data  from  this  period  was  made  available  in  com¬ 
pressed  form  on  magnetic  tapes.  ARCON  was  requested  to  provide  programming 
support  for  the  analysis  of  the  data. 

The  purpose  of  this  effort  was  to  make  recommendations  for  improving 
the  assessment  of  clutter  by  studying  the  clutter  behavior.  A  detailed 
presentation  of  the  background  and  design  of  the  analysis,  with  results  of 
the  initial  phase,  and  proposal  for  further  work,  may  be  found  in  an  in- 
house  report  [Ref.  1  ]. 

The  discussion  below  will  focus  on  the  programming  effort,  in  particular 
on  four  programs.  Each  program  comprises  a  tandem  set  of  subprograms,  the 
first  of  which  performs  the  required  data  sort.  The  remainder  provide 
options  for  presenting  the  data  sort  in  printed  displays. 

1.  POWERLEVEL / POWERD I SPLAY 

2.  BINSORT/BINDISPLAY 

3.  SPEEDSORT/SPEEDDISPLAY 

4.  SPREADSORT/ SPREADDISPLAY 

The  summary  data  is  stored  in  unformatted  (binary)  form  on  magnetic  tapes, 
separately  by  season.  Each  data  record  contains  a  header  with  information  such 
as  date  and  time,  transmitter  power,  operating  frequency,  wave  repetition  fre¬ 
quency,  and  range-azimuth  values.  The  corresponding  range-azimuth  array  fol¬ 
lows,  in  which  the  32  cells  for  each  combination  represent  the  doppler- 
shifted  power  return  to  the  radar. 

Figure  4-1  shows  a  sample  record.  Note  that  the  actual  ERS  array  element 
is  the  negative  of  the  printed  element. 
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4.2  Program  POWERLEVEL 


For  each  spectrum,  the  program  POWERLEVEL  counts  the  occurrence  of  power 
above  a  specified  value  for  each  range-azimuth  element  in  three  doppler  bin 
groups.  The  counts  are  separated  further  by  operating  frequency,  Kp  value, 
and  time  of  day. 

7-9  15-17  23-25 

H - H-H - H-H - -H-H- - h 

t  tiff  •  i i  i  «iii  i 

I  till  111!  till  I 

H - H-H - H-H - H-H - H 

P>-100  P>  -105  P>-100 

The  counts  are  separated  according  to 


Azimuth 

(7  levels) : 

10  to 

80  deg 

Range 

(8  levels) : 

1380 

to  3000 

km 

Bin  group 

(3  levels) : 

7-9, 

15-17, 

23-25 

Kp  value 

(3  levels) : 

0-6, 

7-15, 

16-27 

Frequency 

(3  levels) : 

6-11, 

12-16, 

17-21 

Mhz 

Time 

(3  levels) : 

0-8, 

8-16, 

16-24 

UT 

The  sorted  data  is  saved  in  a  permanent  file.  P0WERD1SPLAY  reads  this 
file  and  prints  range-azimuth  arrays  for  combinations  of  doppler  group,  Kp 
value,  frequency,  and  time  of  day.  Each  printed  element  includes  the  number 
of  counts  made  and  the  number  of  findings  of  power  above  the  specified  level. 
Another  display  option  provides  data  counts  combined  over  all  Kp  values,  all 
frequencies,  etc. 

Figure  4-2  shows  a  diagram. 


4.3  Program  B INSORT 

For  each  spectrum,  the  program  BINSORT  applies  a  filter  to  the  ground 
clutter  level  (bins  15-17).  A  count  is  made  of  the  number  of  occurrences  of 
average  power  greater  than  a  specified  value.  This  count  is  placed  in  the 
category  Looks,  separated  by  azimuth,  range,  and  time  (24  levels,  each  1  hour) 
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If  the  filter  condition  is  satisfied,  then  the  doppler  shifted  power  average 
is  computed  for  bins  7-9  and  bins  23-25,  respectively.  This  sample  is  placed 
in  one  of  the  categories 

1.  Approach 

2.  Recede 

3.  Flat 

(or  counted  as  "other")  based  on  a  comparison  of  the  power  averages. 

Figure  4.3  shows  a  diagram. 

Filter  and  sort  levels,  as  well  as  the  power  difference  are  inputs  to  the 
program. 

The  sorted  data  is  saved  in  a  permanent  disk  file.  A  set  of  DISPLAY 
programs  prints  the  data  counts  in  range-azimuth  arrays  for  each  category 
(Approach,  Recede,  Flat).  Each  array  element  includes  the  number  of  Looks 
and  the  number  of  findings  for  the  category. 

The  range-azimuth  arrays  of  data  counts  may  be  printed  for  1-hour  time 
intervals  or  for  24  hours  combined.  In  addition,  the  arrays  may  be  printed 
for  4-hour  time  groups,  shifted  by  1-hour  increments: 


Time  Group  1 : 

0000-0400  UT 

2: 

0100-0500 

23: 

2200-0200 

24: 

2300-0300 

Another  variation  provides  the  number  of  "hits"  of  Approach,  Recede  and 
Flat  in  percent  relative  to  the  total  number  of  Looks.  The  4-hour  time  groups 
are  used,  as  above,  but  with  no  separation  into  range-azimuth  arrays. 

4.4  Program  SPEEDSORT 

The  program  SPEEDSORT  implements  a  formula  for  associating  each  bin 
with  a  speed  value,  as  follows: 

v  *  -  VMAX  x  A  [m/s] 
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where 


VMAX  =  (WRF/4)  x  (c/f  ) 

op 


(2) 


WRF  =  wave  repetition  frequency  (Hz) 
f  =  operating  frequency  (Hz) 

Q 

c  =  velocity  of  light  =  3  x  10  m/s 

and 

A  =  [Kq  -  2(K-1)]/Kq  (3) 

K  *  number  of  bins  ®  32 
o 

K  ■  bin  number  (1-32) 

so  that  K  -  1,  2, - -  16,  17,  18,  ....  32 

yields  A  -  1,  +.937,  ....  +.0625,  0,  -.0625,  ....  -.937 

[Note:  currently,  bin  1  is  omitted  to  provide  symmetry.] 

This  calculation  is  performed  for  each  data  record.  A  search  is  then 
made  for  speed  values  within  a  specified  interval  [e.g.,  (-125  to  -75)  or 
(75  to  125)].  The  corresponding  bin  numbers,  if  any,  are  retained. 

For  each  spectrum,  a  filter  is  applied  to  the  power  in  bin  17.  A  count 
is  made  of  the  number  of  occurrences  of  power  greater  than  a  specified  value. 

This  count  is  placed  in  the  category  Looks,  separated  by  azimuth,  range,  and 
time  (24  levels,  each  1  hour).  If  the  filter  condition  is  satisfied,  then 
tl.e  doppler-shif ter  power  average  is  computed.  The  bins  are  selected  ac¬ 
cording  to  their  corresponding  speed  value. 

This  sample  is  placed  in  one  of  the  categories 

1.  Approach 

2.  Recede 

3.  Flat 

or  counted  as  "other"  based  on  comparison  of  the  power  averages  (as  in 
B INSORT). 


The  sorted  data  is  saved  in  a  permanent  disk  file.  Display  options  are 
the  same  as  for  BINSORT. 

The  speed  interval  is  an  input  to  the  program  along  with  the  filter 
and  sort  levels,  as  well  as  the  power  difference. 

4.5  Program  SPREADSORT 

The  program  SPREADSORT  calculates  the  difference  between  the  power  in 
bin  17  and  the  average  power  over  bin  groups  corresponding  to  a  specified 
speed  interval.  That  is,  no  filter  is  applied  at  the  ground  clutter  level 
or  the  doppler-shif ted  level.  The  placement  of  spectra  into  categories  is 
determined  by  the  spread  between  these  levels.  The  bin  groups  are  selected 
according  to  the  speed  formula  described  for  SPEEDSORT.  A  count  is  made  for 
occurrences  of  differences  greater  than  a  specified  quantity  (e.g.,  40  dBw) . 
This  sample  is  placed  in  a  category  as  follows: 

1.  Double  -  if  the  condition  is  satisfied  on  both  sides  of  the  spectrum 

2.  Approach  -  if  the  condition  is  satisfied  on  the  left  side  of  the 

spectrum  ** 

3.  Recede  -  if  the  condition  is  satisfied  on  the  right  side  of  the 
spectrum 

Display  options  are  the  same  as  for  BINSORT  and  SPEEDSORT. 


The  speed  interval  and  power  different’  are  inputs  to  the  program. 


'J1  ■.,_,.B  ",l  '_■  *.-  '.v  -  i  -  .  *1.  '  t  » y y ’•■b^  ;.'• .  - 
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Fig.  4-2 
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5.  Meteor  Scatter  Model 


ARCON  was  requested  to  provide  programming  support  to  demonstrate  the 
operation  of  an  experimental  meteor-scatter  model.  An  outline  of  the  basic 
computation  flow  follows: 


I.  Determine,  from  the  geographic  location  of  the  stations,  the  quantities 

f  and  a.  (For  Thule-Sondrestrom,  f  =  600  km,  a  =  6400  km.)  h^  =  80, 

h  -  120. 

u 

II.  Choose  a  direction  of  meteorite  flow, 

0  -  135°,  4>  -  45°. 

III.  Produce  a  series  of  "b-sheets",  one  for  each  value  of  b  (for  testing, 
b  ■  130  to  400,  in  10  km  steps). 

IV.  For  a  given  b,  x  varies  in  the  interval: 

-b  _<  x  b  [  1  km  steps] . 

Calculate  quantities:  y,  z,  u,  v. 

-ABx  +  /(ABx)*  -  CD 

y  " - c - 

A  *  sin8  cos<£ 
ba 

B  -  ^3  —  sin0  sin<j> 

c  *  b5T-^"fT  f 2  sin20  sin2<]>  +  cos2©) 

D  *  (sin20  cos2<]>  +  cos20)  x  2  -  b2  cos20 

/L  2  b2  '  2  2 

z  “/b  -  b2-~F y  "  x 
u  *  x  sin$  -  y  cos<j> 

v  *  -x  cos0  cos$  -  y  cos8  sin$  +  z  sin0 


Calculate  restrictions. 


Tangent  Plane  ®  R 


f 


(9) 


Tangent  Plane  (a)  T 


z  > 


z  > 


/a  2 t1 

f 

A2  -  f2 


(f-y) 

(f+y) 


Lower  Zone  -  Boundary  z  _>  /  (a  *  h.)2  -  (x2  +  y2)  -  >42  -  f2 


Upper  Zone  -  Boundary  z  _>  «/(a  +  h  )2  -  (x2  +  y2)  -  i^2  -  f2 


(10) 

(11) 

02) 


V.  Plot  a  restricted  hot-line  u  vs.  v  fore  each  value  of  b. 

VI.  The  extrema  of  the  restricted  hot-lines  define  an  area  ("cross-section") 
in  u-v  space.  Measure  this  area  in  square  km,  and  tabulate  vs 

VII.  Repeat,  starting  at  II,  using  a  new  meteorite  direction  0, <J>. 

The  programming  effort  has  progressed  through  Step  V.  Printed  output 
was  obtained  for  each  b  value. 

A  plotting  structure  provides  various  options: 

1.  Plot  u,v  corresponding  to  z  values  that  satisfy  the  given  restrictions. 

2.  Plot  all  u,v;  use  special  symbol  to  indicate  those  in  restricted  area. 

3.  Plot  x,y;  use  special  symbol  to  indicate  those  corresponding  to  restricted 

u,v. 

The  TEKTRONIX  Graphics  Terminal  and  its  Hard  Copy  unit  were  used  for 
plot  output.  Examples  of  the  plotted  output  are  shown  in  Fig.  5-1  -  5-3. 
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6.  Direct  and  Multipath  Communications  Data  Analysis 


6. 1  Introduction 

The  task  was  to  provide  software  support  for  reading,  processing,  and 
plotting  communications  data  representing  signal  amplitude  as  a  function  of 
several  variables.  In  particular  the  work  centered  about  the  investigation 
of  the  following  dependencies: 

(a)  amplitude  as  a  function  of  atmosphere  conditions; 

(b)  amplitude  as  a  function  of  the  transmitting  antenna's  elevation 

(c)  amplitude  as  a  function  of  the  time  between  the  initial  signal 

transmission  and  corresponding  signal  reception. 

Data  was  studied  as  a  function  of  two  variables  by  combining  dependence  (a) 
with  dependence  (b)  or  (c) .  The  first  combination  is  known  as  angle-of- 
arrival  data,  or  AOA  data,  and  the  second  combination  as  time-of-arrival 
data,  or  TOA  data. 

The  data  were  invariably  written  on  magnetic  tapes  with  each  tape  con¬ 
taining  approximately  24  hours  of  data,  i.e.,  each  tape  was  associated  with 
data  for  exactly  one  day  of  the  given  year.  Since  for  every  day  of  the 
Gregorian  calendar,  one  can  find  the  Julian  calendar  day,  we  labeled  the 
tapes  using  Julian  dates.  For  example,  tape  BA2256  represents  angle-of- 
arrival  date  for  1982,  day  #256,  i.e.,  September  13.  Letter  B  was  used 
to  indicate  the  data  tape  was  a  backup  copy  of  the  original  tape. 

6.2  Data  Tape  Content  and  Task  Orientation 

Each  tape  contained  approximately  200  records,  usually  from  150  up  to 
220  records.  Each  record  contained  64  signals,  with  each  signal  display 
referred  to  as  a  signal  scan.  The  length  of  the  scan  -  the  number  of  points 
in  the  scan  -  was  different:  64  for  angle-of-arrival ,  and  192  for  time-of- 
arrival  data.  The  time  interval  between  two  sequential  signals  equaled  6 
seconds. 


The  first  main  task  was  to  produce  perspective  on  (pseudo-)  3-dimensional 
plots.  Because  technical  details  such  as  length  of  axes,  tickmarks,  headers, 
labels  and  so  on  were  different  from  AOA  and  TOA  data,  2  sets  of  programs 
were  prepared.  Thus,  program  AOATHREED  produced  3-dimensional  plots  for 
angle-of-arrival  data,  and  program  TOATHREED  produced  3-dimensional  plots 
for  time-of -arrival  data.  These  programs  used  logarithmically  transformed 
data.  We  also  prepared  slightly  modified  programs,  LINAOA  and  LINTOA,  in 
order  to  accommodate  a  request  to  use  linear  data.  Detailed  description  of 
all  control  cards  and  input  cards  is  given  in  Appendix  A. 

The  next  task  was  to  perform  cross-correlation  analysis  between  angle- 
of-arrival  data.  Program  ADACOR  was  constructed  and  consisted  of  two  logi¬ 
cally  independent  parts:  the  second  part  for  making  3-dimensional  plots,  as 
usual,  and  the  first  part  for  performing  cross-correlation  analysis  via  Fast 
Fourier  Transform  technique  (FFT) . 

In  general,  we  used  the  second  scan  from  the  first  record  of  the  cur¬ 
rent  tape  as  a  standard  signal  and  during  the  cross-correlation  process  all 
other  scans  from  the  same  tape  were  compared  with  this  standard  signal. 

(Note:  Of  course,  it  would  be  even  more  natural  to  use  the  first  scan  as  a 
standard  signal,  but  the  first  scan  was  always  defective.)  We  also  had  to 
provide  for  the  case,  when  the  second  scan  was  not  desired  as  a  standard 
signal.  So,  program  CORWRITE  prepared  a  permanent  file,  containing  a  stan¬ 
dard  signal,  and  program  CORREAD  then  used  this  standard  signal  for  correl¬ 
ation  purposes. 


Correlation  Analysis 


The  cross-correlation  function  of  two  sets  of  random  data  describes  the 


general  dependence  of  the  values  of  one  set  of  data  on  the  other.  Consider 
the  pair  of  time  history  records  x(t)  and  y(t).  An  estimate  for  the 
cross-correlation  function  of  the  values  of  x(t)  at  time  t  and  y(t) 


at  time  (t  +  o)  may  be  obtained  by  taking  the  average  product  of  the  two 
values  over  the  observation  time  T.  For  stationary  data,  the  resulting 
average  product  will  approach  an  exact  cross-correlation  function  as  T 
approaches  infinity.  That  is, 


R 

xy 


(a) 


£im  — 
T 

X-KX) 


T 

* 

x(t)y (t+a)dt. 

« 

o 


(1) 


If  R  (o)  =  0,  x(t)  and  y(t)  are  said  to  be  uncorrelated.  For  a  faster 
xy 

implementation  of  this  process,  we  used  the  Fast  Fourier  Transform  Technique 
instead  of  the  direct  (standard)  method. 


6 . 4  Fast  Correlation 

An  infinite-range  Fourier  Transform  of  a  real-valued  or  a  complex-valued 
record  x(t)  is  defined  by  the  complex-valued  quantity 

QO 

X(f)  -  |  x(t)  e"j27rft  dt  (2) 

—00 

Theoretically,  this  transform  X(f)  will  not  exist  for  an  x(t)  which  is  a 
representative  member  of  a  stationary  random  process  when  the  infinite  limits 
are  used.  However,  by  restricting  the  limits  to  a  finite  time  interval  of 
x(t),  say  in  the  range  (o,T),  then  the  finite-range  Fourier  Transform  will 
exist  as  defined  by 

T 

X(f,T)  *  |  x(t)  e“i27Tftdt  (3) 

o 


The  discrete  version  of  this  equation  is 


X(£,T) 


N-l 

hX!  X°  8XP  E"j27Tfnh) 

n=o 
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The  usual  selection  of  discrete  frequency  values  for  the  computation  of 
X(f ,T)  is 


k  =  0,1 . N-l  . 


(5) 


Fast  Fourier  Transform  (FFT)  methods  are  designed  to  compute  the  quantities, 
X^,  and  can  also  be  used  to  compute  the  coefficients  for  a  standard  Fourier 
series. 


Returning  to  cross-correlation  functions,  there  are  two  approaches  to 
the  estimation  of  cross-correlation  functions,  namely,  the  direct  approach 
and  the  FFT  approach. 

The  procedure  for  the  second  case  is  briefly  described  below. 


Vr,h)  ■  iP?  Y.  Xn  yn+r  1  r’0-1 . .  <6) 

n=l 

(If  N  »  m,  it  may  be  more  convenient  to  divid.e  by  N  instead  of  dividing 
by  N-r.) 

The  initial  sample  size  for  both  x(t)  and  y(t)  is  assumed  to  be 

p 

N*2  .  Here,  x(t)  is  the  standard  time  signal  against  which  all  others 
are  to  be  cross-correlated;  y(t)  is  the  time  signal  to  be  cross-correlated. 

Two  separate  FFTs  are  involved  in  these  computations;  one  for  x(t) 
and  one  for  y(t).  We  computed  the  product  of  the  two  transforms,  then  their 
inverse  transform  again  using  the  FFT  procedure.  It  yielded  an  unnormal¬ 
ized  cross-correlation,  and  the  normalization  was  tailored  to  the  initiator's 
needs. 


In  summary,  we  operated  with  3  sets  of  3-dimensional  programs. 
TOATHREED  and  LINTOA  are  for  tapes  with  time-of-arrival  data.  AOATHRE^.D 


and  LINAOA  are  for  tapes  with  angle-of-arrival  data.  AOACOR,  along  with 
CORWRITE  and  CORREAD,  are  for  obtaining  cross-correlograms  for  angle-of- 
arrival  data  only. 

Plots  for  certain  days  showed  unusual  features  and  led  to  the  making 
of  separate  plots  for  each  of  the  scans,  belonging  to  the  area  of  interest. 
This  led  us  to  make  standard  (2-dimensional)  plots.  Again,  3  different 
sets  of  programs  were  made:  time-of-arrival;  angle-of-arrival;  and  cor- 
relograms,  respectively.  "Single  scan"  plots  were  also  produced. 

Using  the  same  system  for  filenames  as  before,  we  created  3  programs: 

AOASCAN  -  for  angle-of-arrival  data. 

TOASCAN  -  for  time-of-arrival  data. 

AOACORSCAN  -  for  correlograms . 

All  programs  are  written  in  FORTRAN  extended  version  4.  for  the  CDC  6600 
computer.  All  plots  were  made  on  the  CALCOMP  off-line  pen  plotter.  During 
the  debugging  process,  the  TEKTRONIX  display  terminal  and  its  hard  copy 
machine  were  used  to  check  the  pictures. 


APPENDIX  A 


Blue  card 


ATTACH, XX, {program  name} 
ID=...,MR=1 

FTN,  ... 

VSN (TAPEl  — 

({tape  number}) 

REQUEST,  TAPEl, 

JH  ,L,  NORING,  NR. 

({ tape  numbers} ) 

REQUEST,  TAPE78, *Q. 

DISPOSE,  TAPE78, *HR. 

ATTACH , X , THREED- 
2WAY , ID*WIK , MR* 1 . 

ATTACH, PEN, EEPLOT. 
or 

ATTACH , TEK , TEKOFFLINE . 


For  "single  scan"  programs: 

CM  =  150K  (not  including  correlation) 

< 

For  3-D  programs: 

CM  =  230K  for  angle-of-arrival  tapes 
CM  =  250K  for  time-of-arrival  tapes 

Main  program  and  subroutines 
(See  list  of  program's  names) 

Compilation 

Tape  with  data 


m 


HI  for  angle-of-arrival  tapes 
HY  for  time-of-arrival  tapes 


Only  for  EEPLOT  plotter. 

Only  for  EEPLOT  plotter. 

Can  be  excluded  for  "single  scan"  pro¬ 
gram,  but  then  #13  has  to  be  changed, 
i.e.  not  contain  X. 

For  EEPLOT  plotter 

For  TEKTRONIX  display 


012  REQUEST, T40140,*PF. 


#12—  ATTACH,  AAA,  PLOTLIB. 


013  LDSET (LIB=X/TEK) 
LDSET(LIB=X/PEN) 
LDSET  (LIB=X/TEK/AAA) 
LDSET (LIB=X/PEN/AAA) 

015  LDSET  (PRESET=ZERO) 

016  LGO. 

016^  CATALOG, T40 140, 
SCANPLOT, ID*  ... 


For  TEKTRONIX  only 


Only  for  "single  scan"  plotting 
See  013  also. 


For  "3-D"  plotting 


For  "single  scan"  plotting 


For  TEKTRONIX  only 


17-  1  or  0 


ISWITCH  (see  main  program) 

1  -  to  make  transformation  of  data 
0  -  to  use  data  as  is  for  plotting 


017-  1  or  0 


018  NEW 


-1  Actually  we  don’t  use  this  card, 
we  only  read  it. 

019  J7J  -  any  integer  number  IISTEP.  See  main  program 

not  0  If  we  want  to  plot  each  scan,  1ISTEP=1 

If  we  want  to  plot  first  scan,  then  skip 
some  of  them,  then  plot  atain,  we  use 
this  parameter. 


#20 

Example: 

256  18 

00 

256 

18 

05 

Start  time  - 

end  time  card 

or 

12 

00 

In  this  form 

it  means  :  we  want  to  start 

from  the  beginning  of  the  tape  and  plot 
during  12  hours,  or  until  end-of-file 
mark. 


021  1  86  060.  Actually  we  don't  use  this  card,  we  only 

read  it,  for  "single  scan"  program  we 
don't  use  it  at  all. 

#22  6  6.  Actual  and  nominal  spacing  of  data 

signals.  (each  6  seconds  we  received 
next  signal.) 

#23  10.  15.  Length  of  X-axis  for  3-D  plotting. 

Usually  10.  for  time-of-arrival  and 
correlation  routines, 

5.  -  for  angle-of-arrival  routines. 

15.  -  angle  of  tilt  for  x-axis 
For  "single  scan"  plotting  we  don't 
use,  only  read  it. 

#24  Angle  of  arrival  Label  card 

13  Sept.  1982 


#25  0 


Number  of  gaps  in  data. 


Example: 

0  25^  256  17  42  51 


Special  "time  of  the  beginning"  for 
time-of-arrival  tapes  ONLY.  To  know 
this  time,  we  read  it  from  corresponding 
angle-of-arrival  tape 


026  6/7/8/9- 

green  card 


pyi 

ATTACH,  TAPE2, 

Only  for  run  program  CORREAD. 

File 

SIGNAL,  ID-. . . , 

SIGNAL  has  to  be  created  before.  We 

MR-1. 

created  file  SIGNAL  during  the 

tion  of  the  program  CORWRITE. 

execu- 

See  below 

REQUEST,  TAPE2,*PF. 

Only  for  run  program  CORWRITE 

#16^ 

CATALOG,  TAPE2, 

SIGNAL,  ID- _ 

Only  for  run  program  CORWRITE. 

Tape 

CC0340  contains  all  programs: 

AOATREED 

TOATHREED 

LINAOA 

LINTOA 

AOACOR 

TOASCAN 

AO  AS  CAN 

AOACOR  SCAN. 

CORWRITE 

CORREAD 

K 

£ 


L' 

t' 

B 


AOATHREED 

LINAOA 

AOACOR 

AOASCAN 

AOACOR  SCAN 

CORWRITE 

CORREAD 

Tape 


TOATHREED 

LINTOA 

TOASCAN 


CC0340, 


TESTLOS 


contains  all  programs 


D.  Nelson,  Perspective  Plotting  of  Two-Dimensional  Arrays:  Computer 
Program  for  a  Digital  Plotter,  Maryland  State  University,  College  Park, 
Maryland,  USA,  1966. 

System/360  Scientific  Subroutine  Package,  IBM,  USA,  1970. 

J.  Bendat,  A.  Piersol.  Random  data:  Analysis  and  Measurement  Procedures 


Wiley- Interscience,  a  Division  of  John  Wiley  &  Sons,  Inc.,  USA,  1971. 


**  ''A 

\  i 

5  MISSION  b 

\  °f  b 

^  Rome  Air  Development  Center  j 

K  RA VC  plans  and  executes  research,  development,  tut  and  v 

*  4  elected  acquisition  pnogaami  in  support  of  Command,  Control  Q, 

St  Communications  and  Intelligence.  1C3 1)  activities .  Technical  £ 
%  and  engineering  support  within  areas  o f  technical  competence  Q 
S  is  provided  to  ESV  Program  deices  ( POs)  and  other  ESD  $ 

J  elements.  The  principal  technical  mission  areas  one  'S 

k  communications,  electromagnetic  guidance  and  control,  sun.-  § 

H  veillance  o ground  and  aerospace  objects,  intelligence  data  v 
»  collection  and  handling,  information  system  technology,  Q 

it  ionospheric  propagation,  solid  state  sciences,  microwave 
%  physics  and  electronic  reliability,  maintainability  and 
s  compatibility. 


X-t 


